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ABSTRACT 

We present the results of rest-fraine, UV slitless spectroscopic observations of a sample of 32 2 ~0.7 
Lyman Break Galaxy (LBG) analogs in the COSMOS field. The spectroscopic search was performed 
with the Solar Blind Channel (SBC) on HST. While we find no direct detections of the Lyman 
Continuum we achieve individual limits (3er) of the observed non-ionizing UV to Lyman continuum 
flux density ratios, /„(1500A )//„(830A ) of 20 to 204 (median of 73.5) and 378.7 for the stack. 
Assuming an intrinsic Lyman Break of 3.4 and an optical depth of Lyman continuum photons along 
the line of sight to the galaxy of 85% we report an upper limit for the relative escape fraction in 
individual galaxies of 0.02 — 0.19 and a stacked 3<r upper limit of 0.01. We find no indication of 
a relative escape fraction near unity as seen in some LBGs at 2 ~ 3. Our UV spectra achieve the 
deepest limits to date at any redshift on the escape fraction in individual sources. The contrast 
between these 2 ~ 0.7 low escape fraction LBG analogs with 2 ~ 3 LBGs suggests that either the 
processes conducive to high f esc are not being selected for in the 2 5 1 samples or the average escape 
fraction is decreasing from 2 ~ 3 to z ~ 1. We discuss possible mechanisms which could affect the 
escape of Lyman continuum photons. 

Subject headings: cosmology: observations — galaxies: evolution — ultraviolet: galaxies 


1. INTRODUCTION 

Reionization brought an end to the cosmic “dark ages,” 
during which the Universe contained a mostly neutral in- 
tergalactic medium (IGM). This transition is thought to 
be triggered by either the radiation from QSOs or from 
galaxies containing large populations of massive stars. 
Recent determinations of the high redshift QSO lumi- 
nosity function suggest that QSO space densities are too 
low to reionize the universe at 2 > 6 (Siana et al. 2008; 
Jiang et al. 2008). The contribution of galaxies to the 
ionizing UV background (Madau et al. 1999) depends 
upon how many of the Lyman Continuum (LyC) photons 
produced by massive young stars that are not absorbed 
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by neutral Hydrogen atoms or dust grains. Direct mea- 
surement of the escape fraction, / e sc, is impossible at 
the epoch of reionization, because intervening absorbers 
make the IGM opaque to LyC photons. Instead, /esc 
must be measured at lower redshifts (2^3.5) in objects 
that are analogous to the galaxies responsible for reion- 
ization. 

Detection of escaping LyC photons has eluded most 
surveys. Observations of Lyman Break Galaxies (LBGs) 
at 2 ~ 3 (Steidel et al. 2001; Shapley et al. 2006; Iwata 
et al. 2009) suggest that in ~10% of starbursts the es- 
cape fraction is quite large, nearing unity. In contrast, 
there are currently no LyC detections locally (2 f 1), de- 
spite tremendous effort (Leitherer et al. 1995; Giallongo 
et al. 1997; Deharveng et al. 2001; Malkan et al. 2003; 
Siana et al. 2007; Cowie et al. 2008, see Siana et al. 2010 
for a review). One explanation suggested by several au- 
thors is the cosmic average escape fraction evolves with 
redshift (Inoue et al. 2006; Siana et al. 2007, 2010). One 
notable difference between high and low redshift studies 
of the escape fraction is the wavelength range used to 
measure the LyC flux. All previous surveys searching for 
intermediate redshift LyC leaking galaxies have utilized 
broad-band photometry which probe ~ 700A while ~ 3 
studies measure the LyC just below the Lyman limit. 
Probing these shorter wavelengths increases the sensitiv- 
ity to the star formation history. Meaning a decrease in 
star formation within the last 10 Myr would significantly 
lower the flux at 700A compared to 900A weakening 
the LyC limits measured from broad-band photometry. 
With this difference in mind and we have undertaken a 
large spectroscopic program with the Hubble Space Tele- 
scope (HST) to study the escape fraction in luminous 
starbursts at 2 ~ 0.7 (GO 11236; PI=Teplitz). 

Using the Solar Blind Channel (SBC) of the Ad- 
vanced Camera for Surveys (ACS), we have obtained 
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far-ultraviolet (FUV) spectra of 32 starbursts chosen to 
be good analogs of LBGs. The observations presented 
in this paper cover rest-frame ~750-1100 A, and pro- 
vide the first spectroscopic search for LyC photons from 
galaxies at 0 < z < 3, as well as the deepest limits on 
the escape fraction for individual galaxies at any red- 
shift. Slitless spectroscopy not only has an advantage 
over broad-band imaging by being able to probe closer 
to the Lyman limit (up to 880A vs. 700A ), but its high 
spatial resolution makes it possible to study any offset 
of LyC radiation compared to UV emission. The large 
sample size provides enough number statistics to make a 
meaningful comparison with the rare LBG detections. 

While no LyC photons are conclusively detected, we 
use these spectra to place strict limits on the UV 
(~1025A) to Lyman continuum (780 — 860A) ratio and 
infer the limit on the escape fraction of ionizing photons. 
In section 2, we describe the selection of LBG-analog tar- 
gets; in section 3, we present the observations and data 
reduction. The method of measuring limits on the escape 
fraction is presented in section 4. In sections 5 and 6, we 
give the results and discuss their implications. A cosmol- 
ogy of SIm = 0.3,11a = 0.70, and H o =70kms -1 Mpc -1 
is assumed throughout. 

2. SAMPLE SELECTION 

Our aim was to choose galaxies at moderate redshift 
which are analogs of the high redshift starbursts pos- 
tulated to be responsible for reionization. The better 
studied z ~ 3 LBG population is conventionally taken as 
accessible analogs to z > 6 galaxies. However, analogs 
at even lower redshift provide the opportunity for a wide 
range of observations, together with the opportunity to 
study the possible evolution of the escape fraction with 
redshift. Hoopes et al. (2007) demonstrated the simi- 
larity between LBGs and the compact subset of local 
UV Luminous Galaxies (UVLG Heckman et al. 2005). 
These sources share many properties with LBGs includ- 
ing luminosity, size, specific star-formation rate, mass, 
and metallicity. The compact UVLG population is quite 
rare locally (about 1 per square degree at z < 0.3), but 
its numbers increase with redshift. They are identified by 
their UV luminosity (> 10 10 L©) and surface brightness 
( I fuv > 10 8 L e kpc -2 ) Hoopes et al. (2007) . 

The wavelength range of the SBC/PR130L configura- 
tion limits the redshift range in which we can measure the 
LyC to 0.6 < z < 0.85 . Furthermore, SBC prism spec- 
tra are slitless, introducing a smearing effect, whereby 
the spectral resolution is degraded proportional to the 
spatial extent of the target. Thus, we were further re- 
stricted our selection to compact galaxies. 

To choose appropriate targets, we needed access to 
deep UV photometry, good redshift estimates, and a 
measure of compactness and surface brightness. These 
three datasets are available in the COSMOS field (Scov- 
ille et al. 2007). We utilized the UV photometry from 
GALEX (Zamojski et al. 2007) and the HST (F814W) 
imaging (Koekemoer et al. 2007) to provide an accurate 
measure of the size and UV luminosity of the source, 
and the photometric redshift estimates of Mobasher et al. 
(2007). From these data, we select targets whose rest- 
frame UV luminosity, Lpuv , and surface brightness, 
I fuv, are the same as those of LBGs. 

We selected non-point sources with large Lpuv and 


I fuv analogous to LBGs (see Figure 1 ). From these 
we discard a small number (~15) which we find by vi- 
sual inspection to have technical problems not obvious in 
the catalog, such as a close proximity to the edge of the 
GALEXfield of view which suffers from artifacts. There 
are 32 sources remaining which meet our criteria. These 
sources are the brightest UV galaxies at z ~ 0.7, but 
their luminosity (Log Luv ~ 10.4 L©) is similar to that 
of typical LBGs at z ~ 3. An ~ L* galaxy at z > 6 
would fall at the faint end of the same range. The UV 
luminosity of the targets corresponds to a star formation 
rate (SFR) of 8-45 M©/yr, without correcting for dust 
extinction. 

The 32 selected galaxies (see figure 2) share many ob- 
servable properties with LBGs. By inference, they are 
the rare, young, strongly star-forming objects at moder- 
ate redshift that have just reached an evolutionary stage 
which was common at z ~ 3. A caveat to note regarding 
the sample selection is that the UV surface brightness 
our targets was derived using the sizes estimated from 
rest-frame optical ACS images, which are typically larger 
than the sizes measured by rest-frame UV images as in 
the case for the Hoopes et al. (2007) criteria. 

3. OBSERVATIONS AND REDUCTIONS 
3.1. Slitless UV Spectroscopy 

The SBC uses a Multi-Anode MicroChannel Array 
(MAMA) that has no read noise and is not sensitive 
to cosmic rays (Boffi 2007). Targets were positioned to 
land on the lower right quadrant of the detector to mini- 
mize the dark current “glow” , which was noted by Teplitz 
et al. (2006) to increase with exposure time. 

The SBC field of view is 31" x35", with a plate-scale of 
0.032' ''pixel- 1 . The SBC PR130L, which covers 1250A 

to 2000 A with a variable dispersion from 1.65 A pixel 1 

at 1250A to 20. 2A pixel -1 at 1800A is used because it 
is not sensitive to the Lya airglow which dominates the 
background in the PR110L prism. The PR130L prism is 
still sensitive to other airglow lines so the observations 
were taken with the SHADOW constraint reducing the 
airglow by a factor of 10. 

The sources have a median GALEX near-UV magni- 
tude of 22.5 (AB), and required 1-4 orbits in the shadow 
position (1500s/orbit; see Table 1 for details). Direct 
imaging (two 480s exposures/orbit) of our targets at 
1 600 A was also acquired in the same visit with the 
F150LP filter to establish a zero-point of the wavelength 
scale. These data were taken over December 2007 - Jan- 
uary 2009 in Cycle 16. An example of the direct SBC 
image and corresponding 2D spectra is shown in Figure 
3 for a typical starburst galaxy in our sample and the 
target which is an AGN-starburst composite (see section 
5.1 for a discussion). 

The flat-fielding and dark subtraction of the SBC data 
was performed by the HST pipeline. Since this work re- 
lies on the measurement of flux below the Lyman limit, 
which is typically on the order of the sky background, sky 
subtraction must be done with care. The mean value for 
a region of sky, typically 300 pixels by 150 pixels above 
and below the target is measured for each prism expo- 
sure. The average sky background in each individual 
exposure is subtracted from each flat fielded prism im- 
age. The spectra were extracted in PyRAF using the aXe 
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slitless spectroscopy reduction package (ver. 1.6 Walsh 
et al. 2006; Kiimmel et al. 2006), specifically designed for 
HST grism and prism data. 

In the extraction, the position of the galaxy on the di- 
rect image is used in combination with the header dither 
parameters to locate the spectrum on the dispersed im- 
age and assigns a wavelength to each pixel in the spec- 
trum. The spectral trace and wavelength solutions are 
defined with respect to a reference position, which is mea- 
sured by running SExtractor (Bertin & Arnouts 1996) on 
a direct image for each orbit. A master catalog was gen- 
erated for each set of prism images taken in a single orbit. 
Using the master catalogs the aXe software generates 
pixel extraction tables (PETs), which contain a spectral 
description for each pixel in the spectrum. Various sizes 
for the extraction window were explored, with an extrac- 
tion width approximately the size of the object maximiz- 
ing the signal to noise of the spectra. The signal to noise 
was further increased (a factor of 1.4) by using aXe’s op- 
timal weighting extraction, which weights the pixels as a 
function of distance from the spectral trace rather than 
giving each pixel the same weight when generating the 
one dimensional spectra. The wavelength solution over 
1300 to 1700A is accurate to within a few angstroms. 
The spectra were flux calibrated using a known sensitiv- 
ity curve for the PR130L prism (Larsen et al. 2006). 

The flux and wavelength calibrated spectra for each 
prism image were co-added weighting each exposure by 
the rms of the background in each individual exposure. 
This was done to account for increasing dark current with 
subsequent exposures. For an individual target the dark 
current was a factor of four higher in the forth orbit com- 
pared to the first. The uncertainty in the spectral flux 
was estimated differently for the wavelengths regions red- 
ward and blue-ward of the Lyman limit. Blue-ward of 
912A where no signal is measured, and we are essentially 
probing the background. The standard deviation of the 
sky flux between the observed wavelengths typically used 
to measure the LyC in our z ~ 0.7 sample was found to 
be comparable to the standard deviation of the flux in the 
same wavelength region of the actual targets. We there- 
fore use the standard deviation of the flux over 780-860A 
rest-frame (the region used to measure LyC flux limits) 
as a conservative estimate of the uncertainty. Red-ward 
of the Lyman break, where flux is detected the errors 
at each wavelength bin in the final ID co-added spectra 
are the ler weighted standard deviation. The final ID 
spectra are presented in 7. 

3.2. Optical Spectroscopy 

Accurate spectroscopic redshifts are required to iden- 
tify the observed wavelengths corresponding to the rest- 
frame Lyman limit and Lyman continuum. Redshifts 
cannot be measured from the FUV spectra themselves 
since the prism probes blue-ward of Ly-a in the rest- 
frame and the resolution is too poor. The photometric 
redshifts, although sufficient for selection purposes, have 
an uncertainty corresponding to ~200A at z~ 0.7. Op- 
tical spectroscopy for 6 of our targets were taken from 
the public zCOSMOS program 12 (Lilly et al. 2009). An 

12 Based on zCOSMOS observations carried out using the Very 
Large Telescope at the ESO Paranal Observatory under Pro- 
gramme ID: LP175.A-0839 


additional 23 spectra were obtained at the Palomar Hale 
5m Telescope using the Double Spectrograph (hereafter 
Double Spec; Oke & Gunn 1982) in February 2008 and 
2009 with seeing conditions ranging from 1” to 2.5” for 
the majority of the observations and light cirrus. This 
instrument uses a pair of CCD cameras which simulta- 
neously obtain long-slit spectra over a “blue” range of 
3500A to 5600A and a “red” range of 6000A to 8500A . 
The observing strategy consisted of acquiring a series of 
600-900s exposures with a slit width of 1.5”. We ob- 
tained a total integration time of 40-90 mins for each of 
the 23 galaxies. 

The data were reduced using standard IRAF tasks 
and wavelength calibrated with He, Ne, and Ar refer- 
ence lamps. The spectral features used for the redshift 
determination were typically a combination of 0[II] at 
3727A , 0[III] at 5007A and in some cases 11,9 at 486lA . 
The spectroscopic redshifts confirmed that all the galax- 
ies were within the 0.6 < z < 0.85 range needed to mea- 
sure the escape fraction with the UV observations. The 
mean spectroscopic redshift of the sample is z=0.679. 
Due to poor observing conditions three targets remain 
with no spectroscopic redshift confirmation. Since only 
one of these (C-UVLG-15) had a detectable UV spec- 
trum and the rest of the sample had photometric red- 
shifts within a few percent of the spectroscopic redshift 
we assumed the photometric redshift for this galaxy. 

4. DERIVING THE ESCAPE FRACTION 

Typically, the escape fraction of Lyman Continuum 
photons is measured relative to the number of photons 
escaping at A reBt =1500A (Steidel et al. 2001). This al- 
lows a straightforward calculation of the escape fraction 
with only two flux measurements (at 1500A and in the 
Lyman Continuum) using the following equation: 


A 


esc,rel 


(fl 50 o/ f Lyc)int 
(/1500/ fL y c)obs X ( S ) 


X e-Xp(T IG M,LyC), 


(1) 


where {fisoo/fLychnt and {fi 5 oo/ h y c)obs are the in- 
trinsic and observed Lyman continuum flux density ra- 
tios. The “LyC” refers to the wavelength at which Ly- 
man continuum is measured (780-860A for this study), S 
is the scaling factor discussed below and TiaM,LyC is the 
IGM optical depth for Lyman continuum photons along 
the line of sight to the galaxy. The intrinsic drop between 
the rest-frame FUV (1500A ) and the LyC (700-900A ) is 
highly uncertain and can not be observed. The assumed 
value of (f isoo / f Lyc)int has varied by study from 3 (Stei- 
del et al. 2001; Shapley et al. 2006) to 6 or 8 (Siana et al. 
2007). 

Since we probe LyC radiation at slightly bluer wave- 
lengths at 830A compared with (Steidel et al. 2001; 
Shapley et al. 2006) we used the SEDs of Bruzual & 
Chariot (2003) and estimated the break amplitude for 
/ 1500//830 to be 3.4 based on / 1500//900 = 3. We assume 
a factor of ~1.2 reduction in fLyc/fioOO for the neutral 
hydrogen opacity in the IGM, modeled in the same man- 
ner as Madau (1995) and Siana et al. (2007). The SBC 
spectra do not measure rest-frame 1500 A for our tar- 
gets. So, to calculate the escaping UV photons at 1500A 
we use the flux within our spectrum between rest-frame 
1000-1050A and apply a scaling factor (S) to estimate the 
flux at 1500A (/i 50 o)- The sources in our sample are by 
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selection, blue objects with low levels of extinction. Fits 
to the photometry suggest an E(B — V) ~ 0.3. Assum- 
ing 0.4 solar metallicity, constant star formation with a 
300 Myr old population the scaling factor to go from the 
/1025 to /isoo is 1.5 in f\. The fluxes derived from the 
spectra are consistent with aperture photometry derived 
using the F150LP images. 

Due to the slitless nature of the spectra there is spec- 
tral “smearing” along the dispersion direction propor- 
tional to the spatial size of the object. A different red 
wavelength cutoff was assumed for each object based on 
the galaxy’s size in the SBC direct image to ensure that 
light red-ward of the Lyman limit did not contaminate 
the LyC flux measurement. The sources have typical 
radii of ~7-20 pixels in F150LP corresponding to a red 
cutoff of range of ~820-880A . Furthermore, the sensitiv- 
ity drops sharply at the blue end and therefore consider 
only the regions at rest wavelengths > 780A . The fi- 
nal spectral region used when estimating the LyC flux 
is between ~780-880A with the red cutoff changing as a 
function of galaxy size (shaded regions in Figure 7). 

The continuum is relatively flat around 850A for galax- 
ies that are actively forming stars (Bruzual & Chariot 
2003), like those in our sample. We therefore integrate 
the observed spectrum over several resolution elements 
to increase the signal-to-noise ratio. The wavelength 
solutions of the SBC, are accurate to a few angstroms 
between observed 1300A and 1700A over the SBC field- 
of-view, and the flux calibrations over this wavelength 
range are accurate to approximately 5% (Larsen 2006). 
The flux below the Lyman limit (/s 3 o) is taken to be 
the average flux between 780-880A (again the red cutoff 
depends on the galaxy size). The uncertainty in fs 3 o is 
derived using the following equation, 


/er,830 : 


>/E(A fl >< AA2) 

A A tot 


(2) 


where A f er is the standard deviation of the flux over 
the assumed LyC region, AA is the size of each pixel 
in angstroms (this value changes as a function of wave- 
length), and A A tot is the total wavelength range be- 
ing averaged. The amount of escaping radiation be- 
low the Lyman limit is typically reported using a rel- 
ative escape fraction (defined earlier in this section) or 
through a UV-to-Lvman continuum flux density ratio 
(/i,(1500A)// 1/ (830A). The former measure however, re- 
quires an assumption for the intrinsic Lyman break which 
is not well constrained. In the next section we present 
both the UV to Lyman continuum flux density ratios and 
the inferred relative escape fractions for completeness. 


5. RESULTS 

We find no direct detections of Far-UV flux in our 
sample of 32 z ~ 0.7 LBG analogs. Nine galaxies (C- 
UVLG-3, 9, 10, 13, 14, 16, 22, 31, 32) were not detected 
in the direct F150LP image or PR130L spectra. The 
non-detections are likely due to in part to their large 
size which resulted in UV surface brightnesses below the 
sensitivity of our observations. Another likely contribut- 
ing factor is that 7/9 of the undetected galaxies had the 
largest extinctions within our sample >0.3. 

The observed flux density ratio, /„1500/ /„830, in the 
individual sources range from 20 to 204 with a median 


of 73.5 (3<x lower limits). In order to convert these ratios 
into a relative escape fraction we apply equation (1) as- 
suming an IGM transmission of 0.85 and a value of 3.4 or 
7 for the intrinsic Lyman break (see Table 1). Our far-UV 
sensitivities give f esc ,rei close to zero, with individual 30- 
upper limits ranging from 0.01 — 0.19. To increase our 
sensitivity further, we stacked the non-detections with 
UV sizes < 0".78 in diameter, which corresponds to a red 
cutoff of ~860A. This red cutoff was chosen to maxi- 
mizing the number of galaxies in the stack while probing 
as closely to the Lyman limit as possible. The stack 
was composed of 18 galaxies placing an upper limit on 
fvlSW/fvSW of 378.7 and a f esc ,rei < 0.01 (see Fig- 
ure 6 . In addition to the global stack we separated 
the sample by morphology, stacking the 8 galaxies which 
appeared in the HST F814W filter to be undergoing a 
merger event and a radius of < ".7 (red cut off of 865A ). 
We find 3cr upper limits for /„1500/ /„830 of 223.2 and 
a fesc,rei < 0.02 for galaxy mergers (see Figure 7 and 
section 6.3 for further discussion). 

5.1. LyC Emission: AGN and Star Formation ? 

Flux below the Lyman limit was detected in C-UVLG- 
17 (hereafter T17), however it is unclear what process is 
responsible for the FUV emission (see Figure 4). T17 is 
detected in the X-ray with a luminosity of ~ 2 x 10 43 ergs 
sec" 1 (2-10KeV), and is a Type 2 AGN based on the op- 
tical spectrum. Although it is highly probable that the 
LyC flux is coming from the AGN it is possible that 
some flux has its origin in young massive star formation. 
The COSMOS 30 photometric band (0.14-8/im) best-fit 
SED is a pure starburst, the IRAC colours do not place 
T17 in the locus of optically bright unobscured AGN and 
is non- variable (over a five year period). The spectrum 
shows no break at the Lyman limit, with ionizing flux 
detected down to A ~ 810 A. Below this wavelength, the 
continuum is strongly suppressed, indicative of a fore- 
ground, high column density absorber at z ~ 0.55. This 
source is interesting as an example of leaking LyC from 
a composite AGN star forming galaxy. 

5.2. Foreground Contamination 

Higher redshift (z ~ 3) studies have detected LyC pho- 
tons in approximately 10% of galaxies (Shapley et al. 
2006; Iwata et al, 2009). Vanzella et al. (2010) estimate 
the at ~ 3 there is a 50% probability that 1/3 of the 
galaxies with directly detected LyC emission are the re- 
sult of foreground contamination by low redshift blue 
galaxies within ~l"of the high redshift galaxy. Figure 9 
highlights an example of foreground contamination found 
in our sample. The NUV GALEX magnitude which is 
used in our sample selection has a PSF of ~ 4"which 
encompasses both galaxies. 

The high resolution HST F814W image (Scoville et al. 
2007) reveals two galaxies separated by ~ 1.2"(<10kpc), 
exhibiting irregular morphologies consistent with a close 
galaxy pair. The southern galaxy (hereafter GalA) has 
two compact nuclei and evidence of long tidal tails. The 
northern galaxy (hereafter GalB) exhibits multiple star 
formation knots, an asymmetric morphology and tidal 
tails. The two primary UV bright knots that make up 
GalB are referred to as GalBl and GalB2 as noted in 
Figure 9. 
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Although GalA was the intended target in our program 
the high resolution SBC direct F150LP image revealed 
that GalB dominated the NUV and FUV flux. The 
Spitzer 24 pm detection is centered on GalA correspond- 
ing to an IR luminosity of L[8_iooopm] = 2.3 x 10 11 Lq 
and a SFR of ~ 40 M Q /yr (Kennicutt 1998). GalA was 
not detected in the FUV (likely due to its high dust 
content), however the UV spectrum of GalB2 (Figure 
10) would imply a relative escape fraction, f re i,esc of 
30% ±6% and a /^(1500A)// y (830A)=14 .3, if it was a at 
the same redshift as GalA. Follow-up spectra with Keck 
DEIMOS revealed that GalB was in fact a low redshift 
interloper at z=0.19. The close (~ 1") projected sepa- 
ration of these galaxies is a good example of foreground 
contamination and could explain some the 2 ~ 3 detec- 
tions of Iwata et al. (2009) particularly the ones with a 
~ l" offset between the LyC emission and assumed LBG 
counterpart. 

6. DISCUSSION 

The deep UV spectra presented here achieve the low- 
est individual limits of the ionizing escape fraction at 
any redshift. We confirm the results of shallower studies 
that suggest low escape fractions in moderate redshift 
starbursts. We now present our findings in the context 
of recent studies and discuss some possible mechanisms 
that could be responsible for the apparent lack of z ~ 1 
starbursts with large escape fractions. 

6.1. Comparison with Recent Work 

Recently, Siana et al. (2010) used HST ACS/SBC 
1500 A imaging and obtain a stacked 3<x limit of 
/„(700)//j/(1500) < 0.02. With the inclusion of previous 
studies (Malkan et al. 2003; Siana et al. 2007) they state 
that no more than 8% of star-forming galaxies at z ~1 
have relative escape fractions greater than 0.50. Cowie 
et al. (2008) stacked the GALEX far-UV (1500 A) fluxes 
of a much larger sample (626 galaxies) and obtain a simi- 
lar 3 <t upper limit of /„(700)//„(1500) < 0.012, confirm- 
ing that starbursts at a ~ 1 have low ionizing emissivi- 
ties. All of these studies have used broad-band UV imag- 
ing to probe the Lyman continuum at ~ 700A whereas 
the detections at z ~ 3 have been obtained through spec- 
troscopy (Steidel et al. 2001; Shapley et al. 2006) and 
narrow-band imaging (Iwata et al. 2008) where the Ly- 
man continuum is sampled just below the Lyman limit 
(880 - 910 A). 

This is the first z ~ 1 spectroscopic study to pro- 
vide Lyman continuum measurements much closer to the 
Lyman limit and is therefore more directly comparable 
to these higher redshift studies. Shapley et al. (2006), 
with a sample of 14 LBGs probe a slightly higher UV 
luminosity and detect two objects with large escape frac- 
tions (with LyC-to-UV ratios significantly above our lim- 
its). One of these objects likely has a lower escape frac- 
tion than initially reported as it was not detected in 
deep narrow-band imaging probing the Lyman contin- 
uum (Iwata et al. 2009). Therefore, only one of 14 ob- 
jects has a significant detection. Similarly, Iwata et al. 
(2009), which probe comparable rest-frame UV luminosi- 
ties as our z < 1 sample, detect large escape fractions in 
~ 10% of LBGs and Lyo emitters at z = 3.1. If the LBG 
analogs in our sample have ionizing properties similar to 


LBGs at z ~ 3, we would expect to detect ~ 3 objects 
with large fractions of escaping LyC radiation, however 
we detect zero. Our individual limits are significantly 
better than the studies at z = 3, so if there are signifi- 
cant numbers of galaxies with lower escape fractions (eg. 
fesc.rd ~ 0.20, rather than unity), we would be able to 
detect them. We see no evidence for this scenario in our 
sample. Comparing our findings, which are in agreement 
with all other z < 1 f esc measurements, with z ~ 3 stud- 
ies (Shapley et al. 2006; Iwata et al. 2009) implies that 
the average escape fraction evolves with redshift, but the 
cause of this evolution remains unknown. We proceed by 
exploring possible mechanisms to explanation that lack 
of large escape fraction galaxies at z ~ 1. 

6.2. Selecting Analogs of High f esc LBGs 

When comparing high and low redshift galaxy samples 
there is always some degree of uncertainty regarding the 
true analog nature of the two populations. As discussed 
in section 2 and shown in Figure 1 we selected a sample 
of LBG analogs sharing many of the same properties as 
the z ~ 3 LBGs observed by Shapley et al. (2006) and 
Iwata et al. (2009). Figure 8, further highlights their 
similarities showings the distribution of reddening, stellar 
mass and rest -frame UV luminosity in these sources is 
similar to the distribution in LBGs. The similarity in 
mass together with the UV-optical colors of the UVLGs 
imply that they may still be undergoing an early, major 
episode of star formation, rather than a small burst on 
top of a hidden older population (also see Barmby et al. 
2004, for the same reason applied to LBGs). 

Ultimately, this sample of LBG analogs shares numer- 
ous similarities to the parent LBG population, but since 
only ~ 10 — 15% of LBGs have been observed with sig- 
nificant LyC detections, perhaps this subclass of LBGs 
have other processes at work allowing or aiding in the 
escape of LyC photons. 

6.3. Galaxy Mergers 

Galaxy mergers offer an intriguing explanation for the 
increased escape fraction seen at z ~ 3. Overzier et al. 
(2008) noted that UVLG typically exhibit faint tidal fea- 
tures suggestive of a merger or recent interaction. They 
therefore propose that the super starbursts in LBGs are 
triggered by gas-rich mergers. As galaxies collide, strong 
gravitational and tidal forces can expel long streams of 
stars, and ignite violent starbursts at rates of a few to 
hundreds of M@/yr (Schweizer 1982; Barton et al. 2000; 
Bridge et al. 2010). During a merger, the tidal fields 
distort the galaxies radially, drawing out galactic mate- 
rial into long tails, plumes and bridges (e.g. Toomre & 
Toomre 1972; Mihos & Hernquist 1996). The HI reser- 
voirs can become disturbed, and the neutral gas pulled 
away from the sources of ionizing radiation producing 
low-column density lines of sight (Hibbard et al. 2000; de 
Mello et al. 2008) through the galaxies, in turn allowing 
the escape of LyC photons. Simulations by Gnedin et al. 
(2008) suggest that the escape fraction in major merg- 
ers can be large (f esc ,rei ~ 30%) compared to non-mergers 
( fesc,rei < 10%) along specific lines of sight. Within 
our sample of 32 galaxies, 11 had morphologies consis- 
tent with merger activity. We independently stacked 8 
of these spectra (removing 3 due to their larger spatial 
extent) and find f esc < 2% (3cr upper limit). 
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If merging is a viable mechanism for clearing pathways 
in the ISM for LyC photons, the orientation of the sys- 
tem along the line of sight is also a likely factor, requiring 
a large sample of UV luminous mergers. Therefore, we 
cannot say whether mergers are an important factor as 
our sample size is at present too small. Currently, there is 
a lack of deep high resolution rest-frame optical imaging 
of the LBGs with larger escape fractions, and the inter- 
pretation of UV morphologies remains problematic (Law 
et al. 2007). Future near-IR observations with HST of 
the z ~ 3 LGB leakers will shred light on this hypothesis. 

As discussed earlier in the section, galaxy mergers are 
capable of clearing pathways, exposing UV bright stars. 
If mergers do facilitate the escape of LyC radiation then 
an evolving merger rate, may be responsible for the ob- 
served evolution in f esc . Numerous observational studies 
and simulations haven shown that the galaxy merger rate 
evolves with redshift, going as ~ (1 + z) 2-3 (Gottlober 
et al. 2001; Conselice et al. 2003; Kartaltepe et al. 2007; 
Hopkins et al. 2007; Bridge et al. 2007; Conselice et al. 
2008; Bridge et al. 2010). The factor of 3-4 increase in 
merger rate between z ~ 1 — 3 as seen observationally, 
would also increase the the number of lines-of-sights and 
range of encounter parameters observed in z ~ 3 galaxy 
mergers by the same factor. This would in turn increase 
the likelihood of detecting LyC at higher redshift. 

6.4. Evolving f esc : Size, Mass, & Star Formation 

The LyC escape fraction is limited by the distribution 
of neutral hydrogen along a line-of-sight and likely de- 
pends on galactic parameters. We now consider what 
galaxy properties could evolve with redshift that reduce 
the efficiency of galactic outflows/chimneys in leaking 
LyC photons from luminous galaxies. 

Typical galaxies (including UV bright galaxies) have 
been shown to be 1.5 to 3 times smaller at z ~ 3 than 
there local counterparts (Trujillo et al. 2006; Papovich 
& Bell 2002; Ferguson et al. 2004). Although our LBG 
analog sample was selected to have similar UV surface 
brightnesses as z ~ 3 LBGs (refer to Figure 1), little is 
known about the true optical sizes of LBGs or the dis- 
tribution of gas. The velocities of galactic winds or out- 
flows have been found to be proportional to the SFR in 
LBGs (Ferrara & Ricotti 2006), therefore LBG and LBG 
analogs, having similar SFR should in principle generate 
outflows with similar velocities (a few hundred km s' 1 ). 
However, smaller galaxies would have higher SFRs per 
unit volume, which can result in more efficient galactic 
winds (Veilleux et al. 2005), more easily clearing path- 
ways or “chimney-like” structures, and in turn allowing 
for higher f esc (Fujita et al. 2003). 

With smaller galaxies, and higher-density starbursts 
comes the potential for a larger fraction of stars born in 
very compact star clusters, including super star clusters 
(SSCs). SSCs can have thousands of young (<50Myr) 
stars within a half-light radius of ~10pc (Veilleux et al. 
2005). These extreme concentrations of hot O and B 
stars can greatly impact the state of the ISM driving 
powerful galactic winds (like those seen in M82), open- 
ing channels for the escaping LyC photons. SSCs have 
been detected in the tidal tails (Charlton et al. 2009) 
and outer regions of galaxies, which could explain the 
spatially offset LyC emission (to the optical emission) 
detected by Iwata et al. (2009) in a few z ~ 3 LBGs. 


There is also some evidence that SSCs found in the local 
group of galaxies have top-heavy IMFs (Harayama et al. 
2008), which enhance the efficiency in clearing lines-of- 
sight, due to the larger outflow velocities generated by 
massive stars. 

Another property which has been shown to possibly 
evolve with redshift is the stellar mass-to-dark matter 
ratio. Recent work by Epinat et al. (2009) and Lethaud 
& Capak (2009), have reported that local massive galax- 
ies have more centrally concentrated dark matter than 
higher redshift galaxies (z ~ 0.7) with comparable stellar 
masses. A higher stellar mass-to-dark matter mass ratio 
at lower redshift would result in a larger gravitational 
potential, and in turn the requirement for larger galactic 
winds to achieve the same “porosity” of the ISM (and 
the probability for ionizing photons to escape their host 
galaxies). Thicker disks due to dense concentrations of 
gas have also been shown in to impede or slow the devel- 
opment of galactic outflows, due to a higher gravitational 
potential (Cooper et al. 2008). 

6.5. Intrinsic LyC-to-UV Ratio: A Top-Heavy IMF? 

An evolving intrinsic LyC-to-UV flux density ratio 
could also be responsible for the observed evolution in 
/ esc . The escape fraction is measured by comparing the 
UV (1500A) flux of a galaxy to flux below the Lyman 
limit. One key assumption that is made involves the 
inherent LyC-to-UV ratio of starbursts. Typically, this 
ratio is considered constant with redshift, but if there 
were an order of magnitude larger production rate of the 
LyC relative to the 1500A, flux at z > 3 than at z < 1, 
the observed change in f esc would be expected. 

Iwata et al. (2009) noted that ~ 3 LBGs in their sample 
had spectral energy distributions intrinsically bluer than 
those expected from population synthesis models, assum- 
ing a standard IMF with moderate dust attenuation. 
They show that an intrinsically bluer spectral energy 
distribution (~ 0.3 mags bluer in NB359-R than those 
of Starburst99), in the absence of QSO activity, can be 
produced with a top-heavy initial mass function (IMF). 
They also suggest that a deviation from the Calzetti et al. 
(2000) dust attenuation law, with less dust absorption 
at 900A, compared to that at 1500A, can come close to 
achieving the observed bluer colors of those LBGs. 

It has been shown that the UV luminosity density at 
z ~ 6 may be insufficient to explain the ionized Universe 
at z > 7 unless the IMF allowed for the production of 
more massive stars (Chary 2008). A top-heavy IMF at 
higher redshift could help explain the larger number of 
LyC detections at z ~ 3 as massive stars produce more 
LyC photons, and stronger supernova driven winds. A 
top-heavy IMF requires regions of potentially low metal- 
licity gas which is consistent with the metallicity and 
dust evolution seen from z ~ 3 to 1 Ferguson et al. 
(2002). 

7. SUMMARY 

Lyman continuum photons produced in massive star- 
bursts likely played an important role in the reionization 
of the Universe. However, their contribution depends 
upon the fraction of ionizing radiation that can escape 
the high column density of HI gas surrounding these star- 
forming galaxies. We have presented HST rest-frame UV 
slitless spectroscopy of 32 z ~0.7 LBG analogs in the 
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COSMOS field to investigate the LyC escape fraction. 
These UV spectra have achieved the deepest limits to 
date on the escape fraction in individual sources at any 
redshift. A summary of our results are as follows. 

We find no direct detections of LyC, consistent with 
other low redshift studies. The individual 3<r upper lim- 
its of / J ,(1500A)// l ,(830A) ratio range from 20 to 204 
(median of 73.5) and 378.7 in the stack of 18 galaxies. 
Assuming an intrinsic Lyman break of 3.4 and an opti- 
cal depth of 85% we report a relative escape fraction in 
individual galaxies of 0.02-0.19 and 0.01 in the stack (3 a 
upper limit). There is no indication of the near unity 
escape fractions found at z ~ 3. The striking contrast 
between the nearly zero escape fractions found in all 32 
z ~ 0.7 LBG analogs with the near unity escape frac- 
tions discovered in 10% of the z ~ 3 LBG population 
strongly argues for an evolving escape fraction. It is un- 
clear however if the lack of near unity escape fraction 
detections at low redshift is due to an evolution in f esc 
itself, or that the number of galaxies with large amounts 
of leaking LyC evolves with redshift. Possible causes 
for a change in the escape fraction with redshift involve 
a top-heavy IMF, larger SFR densities, stellar mass-to- 
dark matter ratios, and/or fraction of SSCs at higher 
redshifts. All these mechanisms enhance the efficiency of 
galaxy winds, increasing the porosity of the ISM, facili- 
tating LyC escape. However, if galaxy mergers aid in the 
escape of LyC radiation then an evolving galaxy merger 
rate could account for the high number of LyC leaking 
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galaxies at z ~ 3. The lack of low redshift galaxies with 
escaping LyC could then be explained by the small num- 
ber of galaxy mergers that have been observed below the 
Lyman limit. 

The escape fraction of UV radiation in z ~ 1 luminous 
starburst galaxies is an important quantity to understand 
since it provides insight into the sources (massive star for- 
mation or QSOs) responsible for reionization. Our study 
has presented a robust measure of the f esc in the low 
redshift Universe and suggests that the escape fraction 
in objects that are analogous to the galaxies thought re- 
sponsible for reionization evolves with redshift. However 
future study is required to isolate the cause of this evo- 
lution. 
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log L fuv [L©] 


Fig. 1. — The FUV luminosity vs. FUV surface brightness for our sample of 32 galaxies (circles and triangles) with 0.65 < z < 0.85 in 
COSMOS. The dash-dotted line indicates the region occupied by z ~ 3 LBGs, while the dashed line marks the region occupied by local 
compact UV luminous galaxies, thought to be analogs for z ~ 3 LBGs (Hoopes et al. 2007). The triangles highlight the galaxies that were 
not detected in the SBC prism observations. The open squares indicate galaxies with no ACS imaging, and the ’x’ notes the target found 
to be an AGN-starburst composite (C-UVLG-17). Note that the angular sizes of our galaxies are derived using rest-frame optical ACS 
images, which are typically larger than the sizes measured by rest-frame UV images as in the case for the Hoopes et al. (2007). 


TABLE 1 

SBC Observations op LBG Analogs in the COSMOS Field 


Object ID 

R.A. 

Decl. 

Exp. Time 

Zspec 

NUV a 

/„830 6 

(fv 1500//p830)o5 s fesc,rel C 

fesc,rel d 


(J2000.0) 

(J2000.0) 

(mins) 







C-UVLG-02 

149.94569 

2.50422 

108 

0.667 

22.40 

2.880e-19 

36.90 

<0.108 

<0.223 

C-UVLG-05 

150.08705 

2.30905 

108 

0.688 

22.60 

1.708e-17 

118.94 

<0.034 

<0.069 

C-UVLG-07 

150.21405 

2.36754 

108 

0.669 

22.69 

9.144e-16 

73.45 

<0.054 

<0.112 

C-UVLG-08 

150.42432 

2.03343 

81 

0.685 

22.35 

7.305e-20 

144.16 

<0.028 

<0.057 

C-UVLG-11 

150.15372 

1.84970 

108 

0.670 

22.39 

6.471e-17 

112.37 

<0.036 

<0.073 

C-UVLG-12 

149.85553 

2.55685 

108 

0.665 

22.50 

1.054e-16 

20.69 

<0.193 

<0.397 

C-UVLG-15 

150.27368 

2.55358 

108 

0.715 

22.60 

3.243e-21 

40.49 

<0.099 

<0.203 

C-UVLG-18 

150.60126 

2.71238 

108 

0.676 

22.53 

1.671e-16 

34.02 

<0.118 

<0.242 

C-UVLG-19 

150.48743 

2.15107 

108 

0.684 

22.41 

4.484e-20 

34.73 

<0.115 

<0.237 

C-UVLG-20 

150.51077 

2.75674 

81 

0.675 

22.09 

5.094e-19 

9312 

<0.043 

<0.088 

C-UVLG-21 

149.62891 

2.19084 

108 

0.701 

22.46 

2.034e-16 

161.10 

<0.025 

<0.051 

C-UVLG-23 

149.88934 

2.73473 

27 

0.709 

21.41 

8.128e-19 

203.56 

<0.019 

<0.040 

C-UVLG-24 

150.17247 

2.63495 

108 

0.752 

22.43 

1.660e-16 

36.05 

<0.111 

<0.228 

C-UVLG-25 

149.62904 

1.66577 

81 

0.671 

22.02 

1.182e-20 

74.98 

<0.053 

<0.109 

C-UVLG-26 

149.60760 

2.16737 

108 

0.676 

22.28 

6.302e-19 

121.33 

<0.033 

<0.068 

C-UVLG-27 

149.69514 

2.74294 

108 

0.652 

22.52 

1.015e-19 

82.99 

<0.048 

<0.099 

C-UVLG-28 

150.45580 

1.65284 

81 

0.657 

22.66 

4.344e-19 

33.84 

<0.118 

<0.243 

C-UVLG-29 

149.46870 

2.58616 

81 

0.704 

21.94 

1.005e-18 

66.43 

<0.060 

<0.124 

C-UVLG-30 

149.51488 

1.89228 

108 

0.694 

22.56 

8.146e-17 

50.50 

<0.079 

<0.163 

C-UVLG-33 

150.76311 

2.80412 

108 

0.693 

22.37 

3.299e-19 

124.97 

<0.032 

<0.065 

C-UVLG-34 

150.69376 

2.84451 

108 

0.738 

22.68 

8.607e-17 

62.97 

<0.063 

<0.130 

C-UVLG-35 

149.69348 

2.59921 

27 

0.678 

21.35 

1.490e-18 

72.07 

<0.055 

<0.114 

Stack (18 Galaxies) 




0.685 

22.14 

378.7 

<0.011 

<0.022 


Stack (8 Galaxy Mergers) 




0.678 

22.50 

223.2 

<0.018 

<0.037 



Note. — The values presents are 3<r limits. Objects 3, 9, 10, 13, 14, 16, 22, 31 and 32 were not detected in the FUV direct image, target 
17 was an AGN 

“NUV magnitude from GALEX public release ver. 3 
4 in units of ergs/s/cm 2 /A 

“assumes and intrinsic flux ratio of (/i 50 o// 83 o)int = 3.4 based in part from (Steidel et al. 2001) 

^assumes (/i50o//s3o)ini = 7 (Siana et al. 2007) 




10 


Bridge et al. 



Fig. 2. — Optical images of the 32 z ~ 0.7 LBG analogs. Thumbnails are F814W HST images except those marked with an asterisk 
(#17, 29,33,34) which are CFHT ground-base i ' . The images are 7" on a side, with the object ID noted in the left corner. 
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Fig. 3. — Example of the SBC F150LP direct UV image (left) and two-dimensional spectra (right) of an AGN (UVLG-17) and a starburst 
galaxy (UVLG-08) in our sample. The Lyman limit is noted in each spectra. Flux below 912A is clearly detected in the AGN, while no 
LyC detection is found in the starburst galaxy. Wavelength increases from the left to right covering ~1220A - 2000A observed. 



Fig. 4. — The individual spectrum of galaxy C-UVLG-17. Below the Lyman limit UV flux is clearly detected, showing that if LyC 
photons were escaping at the ~ 15% level it would be detected by our observations. The source of these photons is likely the Type 2 AGN. 
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Fig. 5. — Deep UV slitless spectra of z ~ 0.7 LBG analogs in the COSMOS field. The spectra have been shifted into the rest frame and 
fluxes are in ergs/cm 2 /a/ A . The horizontal lines (red) represent the average flux over two wavelength ranges, ~780-880A and 1000-1050A 
with 3 cr error bars. 
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Fig. 6. — The stacked spectrum of 18 z ~ 0.7 LBG analogs. The 3a upper limits on the observed UV-to-LyC flux density ratio is 378.7. 
Assuming an intrinsic Lyman break of 3.4 the upper limit on the relative escape fraction is 0.01 (3<r). The Lyman limit is indicated by the 
vertical dashed line. The horizontal lines (red) represent the average LyC and /1025 flux used to derive the escape fraction. Error bars are 
discussed in section 3.1 and are the same as in Figure 7 



Rest Wavelength [A] 

Fig. 7. — The stacked spectrum of 8 z ~ 0.7 LBG analogs with a merger morphology. The 3cr upper limits on the observed UV-to-LyC 
flux density ratio is 223.2. Assuming an intrinsic Lyman break of 3 the relative escape fraction is 0.02 (3cr). The Lyman limit is indicated 
by the vertical dashed line. The horizontal lines (red) represent the average flux over two wavelength ranges. 800-860A and 1000-1050A 
with 3 a error bars. 





Escape Fraction at z ~ 0.7. 


15 





Fig. 8. — Histograms of extinction (left), stellar mass (right) and rest-frame UV luminosity (bottom) for the 32 local LBG analogs. The 
darker shaded regions highlight the 32 galaxies presented in this paper. A sample of z ~ 3 LBGs (clashed lines) in all three panels axe 
taken from Shapley et al. (2005). The dash-dotted vertical lines show the average extinction (Verma et al. 2007), and stellar mass (Yan 
et al. 2006) for z > 5 galaxies and the typical L „ for 2 ~ 6 galaxies (Bouwens et al. 2006). (Bottom) The 7 direct detections of LyC flux in 
LBGs reported in Iwata et al. (2008) are denoted by the lighter shaded histogram, and the lined (angled) histogram shows the 2 ~ 3 LBG 
sample of Shapley et al. (2006). The extinction and stellar mass for the COSMOS sources were estimated in the photometric redshift fitting 
(Mobasher et al. 2007), and the UV luminosity from the GALEX NUV (Zamojski et al. 2007). The histograms show that the galaxies 
presented in this paper have low levels of extinction, similar stellar masses and rest-frame UV luminosities compared to 2 ~ 3 LBGs. 
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Fig. 9. — Postage stamp images left to right are: HST F814W zoomed in followed by F814W, SBC FUV F150LP and the FUV 2D 
spectra. Galaxy A (GalA) and Galaxy B (GalB) are noted in the images. The FUV flux and sites of escaping LyC are concentrated in 
knots (highlighted by circles) in GalB and are referred to in the text as GalBl and GalB2. The merging pair is separated by ~10kpc on 
the sky. Boxes approximate the extraction regions of the UV spectra in the 2D spectral image. Flux below the Lyman limit (assuming a 
z=0.772) is seen to the right of the arrow. 
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Fig. 10. — One-dimensional rest-frame UV spectra of GalB2, a low redshift interloper at z=0.19. 


